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Abstract. We describe an approach to the study of phase transitions in Potts mod- 
els based on an estimate of the complexity of the locus of real zeros of the partition 
function, computed in terms of the classes in the Grothcndieck ring of the affine alge- 
braic varieties defined by the vanishing of the multivariate Tutte polynomial. We give 
completely explicit calculations for the examples of the chains of linked polygons and 
of the graphs obtained by replacing the polygons with their dual graphs. These are 
based on a deletion-contraction formula for the Grothcndieck classes and on generating 
functions for splitting and doubling edges. 
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1. Introduction 



It is well known that the partition function of a Potts model with q spin states on a 
graph G is given by the value at q of the multivariate Tutte polynomial of the graph, 
a famous combinatorial invariant of graphs. The problem of phase transitions in Potts 
models is related to the behavior of the sets of complex zeros and real zeros of these 
polynomials, see for instance [25] , 

In this paper we propose a new approach, based on algebraic geometry, and especially 
on motivic invariants such as classes in the Grothendieck ring of varieties, to study how 
the set of zeros of the partition function of a Potts model changes for a nested family of 
finite graphs that grow in size to approximate an infinite graph. 

We aim at estimating the "topological complexity" of the set of real zeros by computing 
its Euler characteristic with compact support, which is also known to provide an estimate 
of the algorithmic complexity of the real algebraic set. In order to compute this invariant 
and estimate its growth over certain explicit families of graphs, we use the fact that the 
invariant is not just topological but also "motivic" , which means that it defines a ring 
homomorphism from the Grothendieck ring of varieties to the integers. 

Thus, we proceed first to compute explicitly the classes in the Grothendieck ring of the 
varieties defined by the zeros of the Potts models partition functions, following techniques 
recently developed to treat a similar problem arising in perturbative quantum field the- 
ory, for algebraic varieties associated to the parametric form of Feynman integrals. The 
varieties arising in quantum field theory can be viewed, up to a duality, as a limit case of 
the ones arising from the Potts model partition functions. 

We first prove an algebro-geometric inclusion-exclusion formula that relates the classes 
for a given graph to those of the graphs obtained by deleting or contracting one edge in the 
graph, and a more complicated algebro-geometric term, which is the variety defined by 
the intersection of the varieties of the deletion and the contraction. This formula is similar 
to an analogous result proved in pQ for the varieties arising from Feynman diagrams in 
quantum field theory. 

We then show that, when iterating simple operations on graphs, such as splitting an 
edge or doubling an edge, the more complicated term in the algebro-geometric deletion 
contraction formula can be simplified due to cancellations in the Grothendieck ring and the 
resulting operation can be described completely in terms of varieties associated to purely 
combinatorial data on the graph. The corresponding recursive relation leads to explicit 
and remarkably simple generating functions for the classes of graphs obtained from a given 
graph by multiple edge splittings or edge doublings. 

We use the formulae obtained in this way to write explicitly the classes in the Grothen- 
dieck ring for the loci of zeros of the partition function on Potts models over graphs given 
by chains of linked polygons. This class of Potts models was already studied by different 
techniques in the literature (see for instance |22j and references therein). Similarly, we 
compute the Grothendieck classes explicitly for similar chains where the polygons are 
replaced by their dual graphs, the banana graphs. 

For these illustrative examples, we then use the expression for the Grothendieck class 
to compute explicitly the Euler characteristic with compact support, and we show that it 
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grows exponentially as the graphs grow in size, thus estimating the corresponding growth 
in complexity of the set of real zeros of the partition function. 

2. Potts models, multivariate Tutte polynomial and hypersurfaces 

We recall some basic facts and terminology about Potts model that we need in the rest 
of the paper. For a more detailed introduction to partition functions of Potts models and 
their relation to graph combinatorics, we refer the reader to |23| . see also the survey [5]. 

2.1. Multivariate Tutte polynomial. The multivariate Tutte polynomial of a finite 
graph G is defined as follows. Let V = V(G) be the set of vertices of G and E — E(G) 
the set of edges. We do not assume G to be connected. One assigns to each edge e £ E a 
variable t e and then considers the polynomial 

(2.1) Z G (q,t)= ]T q k{G "> J] *«» 

G'CG eefi(G') 

where fc(G') = bo(G') is the number of connected components, and the sum is over all 
subgraphs G' C G that have the same number of vertices V(G') — V(G) of G. Each 
subgraph G' therefore corresponds to a choice of a subset A C E(G) of edges of G, so that 
E(G') = A. The variables t e and the additional variable q are commuting variables. 

A detailed account of the relation between the multivariate Tutte polynomial and the 
physics of Potts models is given in the survey [23]. To briefly recall the main point, in 
the case where q is a positive integer, one considers a g-state model on a graph G, where 
each vertex carries a "spin" that can take q possible values (the case q = 2 recovers the 
usual ±1 states of spin). We let 21 be the set of cardinality q of the possible spin states. A 
state of the system is an assignment of a spin state to each vertex and the energy H of a 
state is the sum over all edges of the graph of a quantity that is equal to zero if the spins 
assigned to the two endpoints of the edge are different and equal to an assigned value — J e 
if they are the same. With the notation t e = e 13,7 " — 1, where (3 is the thermodynamic 
parameter (an inverse temperature up to the Boltzmann constant), one has t e > in the 
ferromagnetic case ( J e > 0) and — 1 < t e < in the antiferromagnetic case — oo < J e < 0. 
The partition function of the system is then the sum over all the possible states of the 
corresponding Boltzmann weight e~^ H , with H the energy of that state. This gives 

(2-2) Z G (q,t)= J2 II ( l + t eSa(v),,( W )), 

<T:V(G)->21 e£E(G) 

where the sum is over all maps of vertices to spin states, v and w are the endpoints 
3(e) = {v, w}, and S is the Kronecker delta. 



It was shown by Fortuin-Kasteleyn [2] that (2.2 1 is the restriction to positive inte- 



ger values of q of a polynomial function in q and that this polynomial is, in fact, the 



multivariate Tutte polynomial (2.1). 



2.2. Deletion— contraction relation. As in the case of the ordinary Tutte polynomial, 



the multivariate Tutte polynomial (2.1) satisfies a deletion-contraction relation. Namely, 
given an edge e € E(G), let G \ e be the graph obtained by deleting the edge e and let 
G/e be the graph obtained by contracting it. One has the following formula. 



Remark 2.1. The polynomial ( |2.1[ ) satisfies 

(2.3) Z G (q,t) - Z G ^ e (q,t) +t e Z G/e (q,t), 

where t consists of the edge variables with t e removed. The deletion-contraction relation 



(2.3) covers all cases, including those where the edge e is a bridge or a looping edge. 
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2.3. The problem of phase transitions. Zeros of the multivariate Tutte polynomial 
are of special interest in relation to the problem of phase transitions of the statistical 
mechanical system described by the Potts model. In fact, the partition function Z G (q,t) 
becomes the normalization factor of the probability distribution on the set of all possible 
states of the system, and zeros of Z G (q, t) would signal the presence of a phase transition 
in the system, for a specific choice of parameters J e and q, and for certain values of the 
inverse temperature /3. 

In the ferromagnetic case, the physical case of interest is where the variables t e G K+. 
For q > 1, the polynomial Z G (q, t) does not have zeros in that domain. The antiferromag- 
netic case with — 1 < t e < is more interesting, and various results on zero-free regions 
are given in |16j . 

Even when there are no zeros in the region of direct physical interest, it is well known 
(see |16j . |23j ) that studying the complex zeros of the polynomials Z G (q,t) can provide 
useful information on phase transitions, not for a single graph G itself, but for a family of 
finite graphs G n , such that Goo = U„G ra determines an infinite graph on which one still 
considers a statistical mechanical system obtained as a thermodynamic limit of the finite 
systems. If loci of complex zeros of the Z G (q,t) can approach the real locus in the limit, 
this will result in the presence of phase transitions for the system on Goo. 

Thus, the geometric problem we concentrate on is to understand and estimate how the 
loci of complex and real zeros, respectively, of the Potts model partition function change 
over certain families of finite graphs {G„} as above. 

Our point of view, in this paper, is to approach this problem from an algebro-geometric 
and motivic point of view, inspired by recent developments on motivic properties of the loci 
of zeros of the closely related Kirchhoff graph polynomials in the setting of perturbative 
quantum field theory, see [I]. [2j. [3]. [8]. [9]. [it]. 

2.4. Potts model hypersurfaces. Studying the zeros of the polynomial Z G (q, t), means 
understanding the geometry of the hypersurface defined by the equation Z G (q,t) = 0. 
Since the polynomial is not homogeneous in its variables, it does not define a projective 
hypersurface (unlike the case of the graph polynomials in quantum field theory), but it 
docs define an affine hypersurface, which we refer to as the Potts model hypersurface. In 
fact, we consider two types of hypersurfaces associated to the Potts model, one where the 
parameter q is treated as a variable q G A along with the other edge variables t ei and one 
where one specializes to a fixed value of q. 

Definition 2.2. Suppose given a finite graph G, with set of vertices V(G) and set of edges 
E(G). Let Z G be the hypersurface in affine space 

A #E(G)+i de fi ned by 

(2.4) Z G := {(q, t) G k* E ^ +1 | Z G (q, t) = 0}. 

For a fixed value of q G A, the hypersurface Z Ga in is given by 

(2.5) Z G . q :={tek* E ^\Z G {q,t)=0}. 

The hypersurface Z G ^ q is therefore a slice of Z G with the hyperplane in A# B ( G ' +1 of 
fixed g-coordinate. In the physically significant cases, one wants to study the complex and 
the real zeros of the hypersurface Z Gtq where q G N is a positive integer corresponding to 
the number of spin-states of the Potts model. 

Definition 2.3. For a finite graph G, the virtual phase transitions of the Potts model 
are the real points Z G (R) of the algebraic variety Z G . For a fixed q, the virtual phase 
transitions are the points of the real locus Z G q (M.) of the variety Z Gq . 

We distinguish here between virtual phase transitions (all real zeros of the polynomial 
Z G (q,t)) and the actual physical phase transitions, which would be constrained by the 
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additional requirement that q £ N and the edge variables are t e > in the ferromagnetic 
case, or —1 < t e < in the antifcrromagnetic case. Thus, for example, in the case of 
a finite graph, even if there are no physical phase transitions in the ferromagnetic case, 
there can still be a non-empty set of virtual phase transitions. 



2.5. The analogy with Quantum Field Theory. In perturbative quantum field the- 
ory, the parametric form of Feynman integrals for massless scalar field theories can be 
expressed as a (possibly divergent) integral of an algebraic differential form on the com- 
plement of an algebraic hypersurface defined by the vanishing of a polynomial associated 
to the graph, the first Kirchhoff polynomial given by 

(2.6) <M*) =£![*«' 

TCGegT 

where t = (t±, . . . ,t n ) are variables assigned to the edges of the graph and T runs over 
maximal spanning forests, that is, subgraphs of G with V(T) = V(G), which are forests 
with bo(T) = bo(G). (Note, the terminology "spanning forest" is used elsewhere for what 
we refer to here as "maximal spanning forests".) 

In the literature on motivic aspects of Feynman integrals [5] , [23] , [55] , it is also common 
to consider the related polynomial 

(2.7) * G (t) = e n *»> 

TCGeeT 

where the sum is, as above, over the maximal spanning forests, but the product is on edges 
in the forest, instead of edges in the complement. 

Definition 2.4. We denote by Xq C A# e ( g > the affine hypersurface defined by the poly- 
nomial (2.6 I and by Xq C P# s ( g ) _1 the corresponding projective hypersurface. 

Similarly, we denote by 3^g C A* b ( g J the affine hypersurface defined by the polynomial 



(2.7) and by 3^g C F^ e ^ 1 the corresponding projective hypersurface. 

Remark 2.5. One obtains ^fc(t) from $g(£) by dividing by IIe6_B(G) * e ana - changing 
variables by the transformation t e l/t e . 

For a planar graph this operation relates the Kirchhoff polynomial of a graph with that 
of a dual graph, see the discussion on the Cremona transformation in [3], 

Remark 2.6. It is also well known (see for instance [IB], [13]) that the graph polynomial 
$g(£) °f (2.7 1 can be recovered from the multivariate Tutte polynomial by the following 
operations: 

(1) Clear an overall factor of q k ^ with k(G) = b (G) the number of connected 
components, that is, consider the normalized Potts partition function 

(2.8) Z G (q,t) =q- k ^Z G (q,t). 

(2) Take the evaluation Z G {q, t)\ q —o. This corresponds to a sum on subgraphs G' with 
k(G') = k(G). 

(3) Of this take then the homogeneous piece with the lowest degree in the t = (t e ) 
variables. This corresponds to the sum over maximal spanning forests, that is, to 
the polynomial ^gW- 



In the context of quantum field theory, Tutte polynomials can also occur, for example 
where one considers scalar field theories on noncommutative spacetimes as in [18) . 
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2.6. The tangent cone. Consider an affine hypersurface X c given by the vanishing 
X = {t E A N | P(t) = 0} of a (possibly non-homogeneous) polynomial P{t) whose leading 
term Pk(t) (the term of lowest order in the variables t — (ti)) is of some degree k > 1. 
Then the tangent cone to X at the origin, TC(X) = TCq(X) is also a hypersurface in A N , 
given by 

(2.9) TC(X) = {teA N \P k (t) = 0}, 

the zero locus of the homogeneous polynomial . This corresponds to the normal cone 
MCq{X) for the subscheme given by the origin C X. The construction known as de- 
formation to the normal cone provides a very useful algebro-geometric replacement of the 
notion of tubular neighborhoods of embedded subvarieties (subschemes) , see [15] . 

In the case of a closed subscheme Y C X, one blows up the locus Y x {0} inside IxP 1 . 
One then considers the complement 

X Y ■= Bl Yx{a} (X x P 1 ) \ Bl Y (X). 

One obtains in this way a fibration Xy — > P 1 , whose general fiber (away from 0) is 
naturally isomorphic to X, while the special fiber over zero is the tangent cone TCy(X). 
This has indeed the effect of deforming X to the tangent cone TCy(X). 
For a finite graph G, we denote by P G the homogeneous polynomial 

(2.10) P G (q, t) = leading term of Z G (q, t), 
in the variables (q, t) £ A# Et - G ^ +1 and by V G the affine variety 

(2.11) V G = {(q, t) G A* E ^ +1 I P G (q, t) = 0}. 

Since the polynomial P G is homogeneous, we can also consider the projective hypersurface 
V G C P# B ( G ). We also consider the affine hypersurfaces 

(2.12) V G , q = {te A# E ^\P G (q,t) = 0} 1 

for fixed q. These are not homogeneous, in general, except in the case q = 0. 



We then have the following rephrasing of Remark 2.6 



Lemma 2.7. The variety V G is the tangent cone of the variety Z G at zero. It has a com- 
ponent given by the hyperplane H = {q = 0} with multiplicity equal to k(G) — bo(G) and 
another component W G , which intersects the hyperplane H along the graph hypersurface 

y G . 



Proof. The first statement follows directly from the definition (2.9) of the tangent cone 



The polynomial P G (q,t) as in (2.10) is the sum of terms of lowest degree in Z G (q,t). 
To see what they parameterize, note that if the subgraph G' determined by a set of edges 
A C E is not a forest, then one or more of the edges may be removed from A without 
affecting the number of connected components, i.e., the exponent of q\ while the degree of 
the product Y[ t e decreases accordingly. 

Further, assume that A is a forest. Then 

k(A) + \A\ = W(G). 

Indeed, this is clear HA = 0; and the left-hand-side does not change if we add an edge 
connecting vertices without closing cycles (\A\ increases by 1, k(A) decreases by 1 for each 
such operation). Therefore, all contributions of forests to Z G (q,t) have degree equal to 
the number of vertices of G, and this is the lowest possible degree. 

Thus, the polynomial P G (q,t) collects the contribution of those terms of Z G (q,t) cor- 
responding to subgraphs that are forests with V(G') — V(G) (spanning forests), 

(2.13) P G (q,t)= Yl 9 MG,) II *«' 

G'CG,6i(G')=0,#y(G')=A r eeE(G') 
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The hypersurface Wq is the locus of zeros of the polynomial Qc(q, t) satisfying Pc(q, t) 
q k ( G )Q G (q t t), where q does not divide Qa(q,t). 

The intersection H n Wq is then given by the locus of zeros of the polynomial 

(2.14) Qa(0,t)= II *«> 

G'CG, max forest eeE(G') 



which is the polynomial $g(0 of (2.7). □ 



2.7. The Grothendieck ring of varieties and mixed Hodge structures. The alge- 
braic varieties Xq and J^G associated to Feynman graphs have been studied extensively in 
recent years in terms of their classes in the Grothendieck ring of varieties, see for instance 
0,0, [H], . We will be applying here analogous techniques to the Potts model hy- 
persurfaces Zq and Zq^ and apply the results to the problem of phase transitions. Thus, 
we recall here a few things about the Grothedieck ring of varieties, and a result of p], 
whose analog for Potts models hypersurfaces we prove in this paper and will be the basis 
of our motivic approach to phase transitions. 

The Grothendieck ring Kq(V%) of varieties over a field IK is generated by isomorphism 
classes [X] of smooth (quasi)projective varieties with the inclusion-exclusion relation 

(2.15) [X] = [Y] + [X\Y] 

for any closed embedding of a subvariety Y C X, and with the product structure given 
by[XxY] = [X][Y]. 

In the following, we will be interested in considering the Potts model hypersurfaces 
as varieties defined over C, but we will also be focusing on their real zeros, hence think- 
ing of them as varieties over K. Thus, in the following we simply write K (V) for the 
Grothendieck ring, whenever the arguments do not depend on what field we work over, 
and we will explicitly mention C or R when needed. 

The class [X] in the Grothendieck ring is a universal Euler characteristic for algebraic 
varieties (see [7]), in the sense that any invariant of isomorphism classes of algebraic 
varieties that satisfies the inclusion-exclusion relation and is multiplicative on products 
factors through the Grothendieck ring. These invariants are sometimes called motivic. 

In particular, in the case of complex algebraic varieties and of classes in Kq(Vc), among 
these motivic invariants that factor through the Grothendieck ring we have the usual 
topological Euler characteristic, but also the virtual Hodge polynomials. These will be 
useful to us in the following so we recall briefly the definition. 

The virtual Hodge polynomial of an algebraic variety is defined as 

d 2d 

(2.16) e(I)(x,j/)=^^(W, with e™(X) = £(-l)*ft™(fr c fc (J0), 

p,q=0 fe=0 

where, for each pair of integers (p,q), the term h p ' q (H* (X)) is the Hodge number of the 
mixed Hodge structure on the cohomology with compact support of X. If X is smooth 
projective, then the virtual Hodge polynomial reduces to the Poincare polynomial, with 
e p ' q (X) = (—l) p+q h p ' q (X) being the classical pure Hodge numbers. 

The fact that the virtual Hodge polynomial factors through the Grothendieck ring 
Kq(Vc) of varieties means that an explicit formula for the class of a variety in the 
Grothendieck ring can be used to compute the virtual Hodge polynomial and obtain some 
explicit information on the Hodge numbers and the mixed Hodge structure. 

2.8. Virtual Betti numbers of real algebraic varieties. As we mentioned above, in 
the case of complex algebraic varieties, the ordinary topological Euler characteristic is a 
motivic invariant. This is not true for real algebraic varieties, as the additive property 
over closed embeddings need not be satisfied. However, it is known that there is a unique 
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motivic invariant that agrees with the topological Euler characteristic on compact smooth 
real algebraic varieties and is homeomorphism invariant (but not homotopy invariant), see 
[2"T] and also [T3], It is defined, for any real (semi) algebraic set S, as 

(2-17) Xo(S)=^2(-l) k b^ M (S), 

k 

where b^ M (S) arc the Borcl Moore Bctti numbers, namely the ranks of the relative ho- 
mologies Hk{S, oo), where S is the Alexandrov compactification. Equivalently, they are 
the ranks of the cohomology with compact support H*(S). 

Example 2.8. Let L = [A 1 ] be the Lefschetz motive, the class of the affine line in Kq(V) 
and let T = [G m ] = L — 1 be the class of the multiplicative group G m = A 1 \ {0}. The 
topological Euler characteristic x '■ Kq(Vc) ~ > Z satisfies x(^) — 1 an d xCO = 0j while 
the Euler characteristic with compact support Xc ■ -Ko(Vr) — > Z satisfies Xc(L) = — 1 and 
Xc(T) = -2. 

Moreover, it is shown in that the Betti numbers with Z/2Z coefficients bk(X) = 
dim Hk{X, Z/2Z), defined in the usual way for compact smooth real algebraic varieties, 
extend in a unique way to Xo(Vk), so that one obtains a ring homomorphism 

(2.18) : K {Vi) ^ Z[t], 

such that 0(X, t) = J2k^k(X)t k for compact smooth varieties. The coefficients 0k of 
the ring homomorphism are called the virtual Betti numbers. They are not topological 



invariants. However, they compute the Euler characteristic (2.171, namely, 
(2-19) 0(X,-l)=Xc(X) 

for all real algebraic varieties X, with both equal to the ordinary Euler characteristic x(X) 
in the compact smooth case. Notice that, while Xc(X) is the alternating sum of the ranks 
of the Borel-Moore homologies, the virtual Betti numbers 0k(X) are in general not equal 
to the Borel-Moore Betti numbers b^ M (X) (for instance, the 0k(X) can be negative), 
although their alternating sums agree. 

In the case of a compact smooth real algebraic variety, which is the real locus X(R) 
of a smooth projective complex algebraic variety X(C), there are ways to bound the 
"topological complexity" of X(K) in terms of invariants of X(C), in the form of Petrovskii- 
Oleifiik inequalities: for example, for X(IR) a smooth projective real algebraic variety of 
even dimension n = 2p, one has |3], [17] 

(2.20) \x(X{R)) - 1| < h^ p (H n (X(C))). 

This type of result was extended to cases with isolated singularites in [TS] , where one gets 
(2.21) 

Eo<,< p ft M (^(^(C))) n = 2p 

. Eo< 9 <^ 9 ' 9 (^ l (^(C))) + ^ +1 ' p+1 (ff c p+1 (^(C)) n = 2p + l. 

However, more generally one does not have a Petrovskii-Oleinik type inequalities to 
estimate the virtual Betti numbers and the Euler characteristic Xc(X) of arbitrary real 
algebraic varieties in terms of the virtual Hodge polynomials of the complex variety. For 
complex varieties the virtual Bctti numbers can be computed in terms of the virtual Hodge 
polynomial. In fact, one can introduce the weight fc Euler characteristic given by setting 

w^(X(C))= ]T h^{H k c {X{C))), 

p+q=j 

which equal b^{X) for j = k and zero otherwise in the smooth projective case, and 
are otherwise equal to the ranks of the quotients of the weight filtration ^(^(C)) = 



\ X (X(R)) - 1| < 
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dim c W/(X)/WjLipf(<C)) on the cohomology with compact support H*(X(C)), Then 
for arbitrary complex algebraic varieties the virtual Betti numbers are given by (20J) 

p j (x(c)) = (-iy'£(-i) k w*(x(c)). 

k 

In general one does not have a good way to estimate the virtual Betti numbers /3fe(X(R)) 
of a real algebraic variety, nor their alternating sum % C (X(R)), in terms of the virtual 
Betti numbers j3 k (X(C)). 

Although in general one cannot estimate x c (X(R)) in terms of a suitable Petrovskii- 
Olemik type inequality, we will show that in certain cases one can compute explicitly 
both Xc{X(M)) and the virtual Hodge numbers of X(C) as a consequence of being able to 
compute explicitly the class [X] in the Grothendieck ring of varieties. 

We will discuss later how these considerations relate to the problem of phase transitions 
in Potts models. In particular, we will see that, by working with classes in the Grothendieck 
ring, we obtain some estimates on the topological complexity of the set of virtual phase 
transitions of the Potts model over certain families of Unite graphs G n approximating some 
infinite graph G — U n G n . 

2.9. Grothendieck classes of Potts model hypersurfaces. We also introduce the fol- 
lowing notation for the classes in the Grothendieck ring of the Potts model hypersurfaces. 

Definition 2.9. Let [Zq] be the class in the Grothendieck ring K Q (V) of the Potts model 



hypersurface (2.4). Also let {Zq} denote the class of the hypersurface complement, 

(2.22) {Z G } = [A* E ^ +1 \ Z G ] = h* E ^ +1 - [Z G ], 

where L = [A 1 ] is the Lefshetz motive (the class of the affine line). The classes [Zo, q ] 
{Zc q } = [A# E ( G > \ ZG,q] are similarly defined for the hypersurface Zc. q of (|2.5[). 



As in the case of the graph hypersurfaces of Feynman graphs (see pQ, [2], 0), we will see 
that for Potts models it is simpler to write explicit formulae for the class of hypersurface 
complement {Zq} than for the class [Zq] of the hypersurface itself, though the information 



is clearly equivalent due to the simple relation (2.22) between them 



2.10. The Grothendieck class for fixed q and the fibration condition. We discuss 
here the relation between the classes of the hypersurface complement {Zq} and {Zq^} 
for the full Potts model hypersurface and for the one with fixed q. We identify a useful 
condition, according to which the the class of {Zq} behaves as one would expect in the 
case of a fibration on the locus q ^ 0, 1. We will later identify specific classes of graphs we 
want to work with and check that they satisfy this condition. We do not address in this 
paper the question of how general this condition actually is, nor the question of whether 
the variety Zq itself really is a locally trivial fibration over the locus q ^ 0, 1, at least for 
some specific families of graphs. 

One can see directly from the poynomial Zc(q, t) why q = and q = 1 should certainly 
be special values, for the following reasons. Recall that the equation we are dealing with 
is 

Z G (q,t)= £ «*< A >n*" 

ACE(G) aEA 

• For q = and G nonempty, this is 0: indeed, k(A) > for every subset A of edges 
in that case. This just says that the hypersurface Zq{<jC) = has a component 
along q — 0; this component can be removed (dividing Zc(q,t) by q k ^) and 
the residual hypersurface may be studied over q = 0. This is the hypersurface 



Qc(0,t) — of (2.14), which is the dual 3^g °f the graph hypersurface Xq as in 



Definition 2.4 This falls back on the case investigated in [2], pQ. 
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• For q = 1, the polynomial becomes 

(2.23) Z a (l,t)= £ II> = 11(1 + *e)- 

ACB a£A e£E 

This is a union of normal crossing divisors, and in fact it consists of essentially n 
coordinate hyperplanes in A". The complement is the set of n-tuples (ti,. . . ,t n ) 
(n = #E(G)) with each U + 1 ^ 0, a copy of the n-torus. Thus the class of its 
complement is T™. 

The condition that the class {Z G } behaves as in the case of a fibration over the locus 
q ^ 0, 1 can then be formulated as the condition that 

• The class {Zc, q } is independent of q for q ^ 0,1; this class will be denoted 

{-Zg, ? ^o,i}; 

• The following holds: 

(2.24) {Z G } = (T — \){Z GtM } + T* E ^ G \ 

This accounts for the fact that the complement of Z G = is contained in q ^ 0, has 
a torus fiber over q = 1, and (heuristically) fibers over q ^ 0, 1, a copy of T — 1, with 
constant fiber class. 

In particular, a necessary condition for (2.24) is that (T — 1) divides {Zq} — T# B ( G ' in 



the Grothendieck ring, so any counterexample to this property would give examples where 



the fibration condition (2.24) does not hold 



3. Deletion-contraction for classes in the Grothendieck ring 

In p] it was shown that, in the case of the graph hypersurface complements A^ E ^\X G , 
the classes in the Grothendieck ring satisfy an algebro geometric analog of a deletion- 
contraction relation. More precisely, it was proved in pQ that, for a graph G with n = 
#E(G) and for a given edge e € E(G), the classes of the varieties X G , X G / e and X G ^ e are 
related by 

(3.1) [A" \ Xq] = L [A"" 1 x (Afe Na n X G/e )] - [A"" 1 \ Af Gxe ], 

when e is neither a bridge nor a looping edge, and 

[A" \ Xq] = L [A"" 1 \ X G/e ] = L [A"" 1 \ X G ^ e ] 
when e is a bridge and 

[A" \ X G ] = (L - 1)[A™ _1 \ X G/e ] = (L — l)[A n_1 \ X G ^ e ] 
when e is a looping edge, where L = [A 1 ] is the Lefschetz motive, as above. 



Notice how (3.1) is not a combinatorial deletion-contraction formula: indeed, the term 
involving the intersection X G ^ e n X G / e of the hypersurfaces of the deletion and the con- 
traction is in general difficult to control explicitly, even if one has an explicit formula for 
the classes of the deletion and the contraction separately. However, it was also shown in pQ 
that, for certain families of graphs, such as chains of polygons, one obtains recursive rela- 
tions in the Grothendieck ring, where the "problematic" term in the deletion-contraction 
formula cancels out and one obtains an explicit generating function for the classes of the 
varieties associated to the family of graphs. The result provides a way to control how the 
class in the Grothendieck ring grows in complexity when the graph is enlarged through 
some simple operations, such as doubling edges or splitting edges. In the setting of quan- 
tum field theory the families of graphs obtained through such simple operations, like the 
chains of polygons, are typically not complex enough to capture interesting behaviors of 
the associated periods, but we will argue here that analogous operations performed in the 
setting of Potts models already gives rise to interesting non-trivial cases. 
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3.1. Algebro-geometric deletion contraction for Potts model hypersurfaces. 



We prove here an analog of the deletion-contraction formula (3.1) for the classes in the 
Grothendieck ring of the Potts model hypersurfaces. We first analyze the case of the full 
Z G and then the case of Z G g with fixed q and of the tangent cone V G at zero. 

We now state our main result on the deletion-contraction properties for Potts model 
hypersurfaces. Notice that, in the following statement, there is no distinction between 
the case of bridges or looping edges and all the other edges, just as in the combinatorial 
deletion-contraction relation for the multivariate Tutte polynomials. 



Theorem 3.1. Let G be a finite graph and e an edge of G. Then the class {Zq} of (2.22) 
satisfies 

(3.2) {Z G } = L{Z G/e n Z G ^ e } - {Z G/e }. 

Proof. The result follows from the combinatorial deletion-contraction relation for the mul- 
tivariate Tutte polynomials 

Z G (q, t) = Z Gse (q, t^ ) + t e Z G/e (q, £« ) , 

where is the set of variables t — (t e ') e 'eE(G) with the variable t e omitted. We then 
check the various cases. 

• If Z G / e (q,t^) 7^ 0, then Z G (q,t) is guaranteed to be ^ if t e does not equal 
—Z G ^ e (q,i^)/Z G / e (q,i^). This accounts for a G m worth of t e 's for each such (q,i^ e '), 
contributing a class 

(L - l){Z G/e }. 

• If Z G/e (q,i^) = 0, then Z G (q,t) ^ if and only if Z G ^ e (q,i^) ± 0. This accounts 
for an A 1 worth of t e 's for each (g,< (e) ) such that Z G/e (q,i {e) ) = and Z G ^ e (q,i^) ^ 0. 
This contributes a class 

l- [z G/e \ (z G/e nz G ^ e )}. 

Note that 

[z G/e ] - [Za/e n z G ^ e ] = h# E ^ - [z G/e n z G ^ e ] - h* E ^ + [z G/e ] 
= {z G / e n z G ^ e } - {z G / e }. 

Thus, the two contributions add up to 

{z G } = (L — i){z G/e + l ({z G/e n z G ^ e } ~ {z G/e }) = h{z G/e n z G ^ e } - {z G/e }, 

as claimed. □ 

The following properties of the classes of Potts model hypersurfaces are simple conse- 
quences of the definitions, or follow easily from Theorem |3.1| 



Corollary 3.2. The classes {Z G } £ Kq(V) of (2.22) satisfy the following properties: 

(1) If G consists of a single vertex and no edges, then {Z G } — L — 1. 

(2) IfG consists of a single edge, with either one or two vertices, then {Z G } = (L — l) 2 . 

(3) If a graph G' is the union G' = G\ U„ Gi of two graphs joined at a vertex v, and 
G" denotes the disjoint union of the same two graphs, then {Z G >} = {Z G n}. 

(4) If G is obtained by joining G\ and Gi with a single edge from a vertex of G\ to a 
vertex 0/G2, then {Z G } = (L — 1){Z G '}, with G' — G\ G2 as above. 

(5) If G is obtained from a graph G by appending a single (looping or otherwise) edge 
to a vertex, then {Zq} = (L — 1){Z G }. 
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Proof. (1) If G consists of a single vertex, then Z G = q defines a point in the affinc line 
A 1 with coordinate q. 

(2) If G is a single edge joining two distinct vertices, then Zc(q 1 t) = qt + q 2 and 
therefore {Z G } = L 2 — 2L + 1 = (L - l) 2 . 

If G is a single looping edge, then Z G (q, t) = qt + q and again {Z G } = (L — l) 2 . 

(3) If G' consists of two graphs G\ and G 2 joined at a vertex, then Zqi (q, t) — ~Z Gt Zq 2 . 
If G" consists of the disjoint union of G\ and G2, then Zo"(q,t) = Z Gl Z Gl . It follows (if 
none of the graphs is empty) that {Zqi} — {Z G n}, and that {Zg> n Z G u} = {Z G >}. 

(4) Applying Theorem |3.1| this implies that if G is obtained by joining G\ and G2 with 
a single edge from a vertex of Gi to a vertex of G2, then {Zq} = (L — 1){Z G i}. 

(5) Let G' be the disjoint union of G and of a single vertex. Then Zqi — qZ G , and 
in particular {Z G i} = {Zq}, and Zqi n Zq = Zq. Applying Theorem |3.1| we see that 
if G is obtained by appending a single non-looping edge to a vertex of G, then {Z^} = 
(L — 1){Z G }. The same conclusion is reached if G is obtained from G by adding a looping 
edge. □ 

Notice that, unlike what happens with the graph hypersurfaces of Feynman diagrams 
(see [1]), in the case of a disjoint union of graphs the class of {Zq iU q 2 } is not the product of 
the classes of the hypersurface complements of the two graphs, since here the polynomials 
Z Gl and Zq 2 have the same variable q in common. However, if the graph G = G\ U G2 is 
given by the disjoint union of two graphs G\ and G2, and all the graphs involved satisfy the 



fibration condition (2.24), then we can obtain an explicit formula for the class {Zq iU q 2 }. 



Corollary 3.3. Let G — G\ U G2 be the disjoint union of two finite graphs G\ and G2, 



and assume that the classes {Zq} and {Zq^ satisfy the fibration condition (2.24 1. Then 
the class {Zq} can be expressed explicitly in terms of the classes {Zq ± } and {Zq 2 } by 
(3.3) 

, {Z Gl }-{Z G2 }-T# E ^ -{Z G2 }~T* E ( G ^ . {Z Gl } + T* e ^)+# e ^ +1 
{2q iU g 2 } - J^l ' 

Proof. For q fixed, the remaining variables are indeed distinct for disjoint G%, G 2 . Thus, 
the classes satisfy 

(3-4) {Z GlUG2 J = {Z Guq }-{Z G2 J. 

(This holds for all q, including the special values q = and q = 1.) If the classes {Z Guq } 



are independent of q ^ 0, 1 and the formula (2.24 1 holds, then we get 

{Z GlUG2 } = (T — l){Z GlUCWOil } + T #E(G l) + # E { G 2) 

= (T — \){{Z GlM } ■ {Z G2 , M })+T* E(G ^ + * E ^ 

= (T 1) ( ^}- T#g(Gl) . {Zg 2 }-T* e ^ \ t#e{Gi)+#e{G2) 
\ ) y T — 1 T - 1 J 

_ {Z Gl } ■ {Zq 2 } - T# E < G ^ ■ {Zq 2 } - T# E ^ ■ {Z Gl } + T #E(G l ) + #E(G 2 ) + l 

T - 1 

□ 

3.2. Deletion— contraction for fixed q. Deletion-contraction works exactly as in the 
case of the full Potts model hypersurface Zq. 

Proposition 3.4. For a finite graph G, the class {Z Gtq } for fixed q satisfies 
(3-5) {Z Ga } = (T + l){Z G/e<q n Z G ^ q } - {Z G/e J. 



The argument is identical to the one used in the proof of Theorem 3.1 (For q = 0, all 



classes equal 0.) We then also have the analog of Corollary 3.2 
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Corollary 3.5. For q 0, the classes {Z Gq \ satisfy the following properties: 

(1) For G a single vertex, {Z Gq } = 1. 

(2) For G a single edge joining either one or two vertices, {Zc, q } = T. 

(3) If G' consists of two graphs G\ and G2 joined at a vertex, or disjoint, then 

{ZG',q} = {ZGuq} ■ {Zo 2 ,q}- 

(4) // G is obtained by joining G\ and Gi with a single edge from a vertex of G\ to a 
vertex of G2, then {Z Gq } = T{Z G > q }. 

(5) Appending a single edge to a vertex of G multiplies the class {Z G ^ q } by T. 

Proof (1) For G a single vertex, Z G (q,t) = q ^ 0; thus Z Ga = C A , and {Z Ga } = 1. 

(2) For G a single edge joining two distinct vertices, Z G (q, t) = qt + q 2 defines a point in 
A 1 ; thus {Z Gq \ = T. For G a single looping edge, Z G (q, t) = qt + q and again {Z G ^ q } = T. 

(3) This case is simpler than the case for {Z G } with indeterminate q. It follows, as in 



( 3.4 ) , from the fact that the polynomials Z Gl (q, t) and Z G2 (q, t) have none of the variables 
other than the fixed q in common. 

(4) This again follows from the simpler formula for unions as in (3), and from the 
computation for a single edge. The relation is the same as in the indeterminate-g case. 

(5) This is again the same formula as in the free case, which here follows from (3) and 
the case of a single edge. □ 

3.3. Deletion— contraction for the tangent cone. In the case of the tangent cone V G 
at zero of the variety Z G , which, as we have seen interpolates between the Potts model and 
the graph hypersurfaces considered in the quantum field theory context, it is convenient 
to introduce the following notation for the classes in the Grothendieck ring. 

We still denote by [V G ] the class of V G and by {V G } = [A# E ^ +1 \ V G ] the class of 
the complement. We also use the notation y G for the graph hypersurface given by the 
intersection of the component W G = {Q G {q,t) = 0} with the hyperplane H — {q = 0}. 
Thus, y G is the locus of zeros of Q G (0,t), where P G (q,t) — q k ^Q G (q,t). This gives, at 
the level of the classes 

(3.6) [V G ] = [W G ]+L #£(G) -[3H 

where L# E( - G ^ = [H] and V G = W G U H with y G = W G D H. Thus, we obtain the 
following. 

Lemma 3.6. The class {V G } of the complement of the tangent cone us 
given by the class {V G } = {W G } — {3^ G }- 



Proof. Using (3.6 1 and taking complements, we obtain 



L #£(G) + 1 _ [Vg] = L #S(G) + 1 _ [Wg] _ h #E( G ) + { y a] = {Wg} _ { y c} _ 

Alternately, one may simply observe that the complement of W G is the disjoint union of 
V G and cY G . □ 

We can then see directly that the class {V G } satisfies the following simple properties. 

Lemma 3.7. The class {V G } satisfies: 

(1) If G' is obtained by attaching a looping edge to a vertex of G, then 

{V G ,} = (Y + l){V G }. 

(2) // G' is obtained by appending a non-looping edge to a vertex of G, then 

{v g ,} = t{v g }. 

(3) More generally, if G is obtained by connecting two disjoint graphs by a bridge e, 
then{V G }=T{V G/e }. 
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(4) If G is obtained by attaching an edge parallel to one of the edges of G, then 

OM = (t + i){v g }. 

Proof. (1) The polynomial P G is only counting forests and loops are excluded from forests, 
so the polynomial P G i equals P G , but it is viewed in one dimension higher. This simply 
multiplies everything by L = T + 1 . 

(2) Attaching an unconnected edge e multiplies P G by (q+t e ). The condition (t+q)P G 

implies Pq ^ and t + q^O. This says that the complement to P G > — fibers over the 
complement to Pq = 0, with T fibers. 

(3) The same argument proves this assertion. 

(4) Let e be the edge of G that we are doubling, and call / the new parallel edge. Then 
Pg> is obtained from Pq by replacing t e by t e + tf. This operation amounts to taking a 
cylinder, hence multiplying everything by L = T + 1. □ 

We then look at the deletion contraction formula. In the cases of {Z G } and {Z G ^ q } 
considered above, we did not have to make a special case for bridges and looping edges 
because the combinatorial deletion-contraction formula for the multivariate Tutte poly- 
nomial does not make such a distinction. However, in the case of the tangent cone, as in 
the case of the graph hypersurfaces of Feynman diagrams of jT] , we need to take these two 
special cases into account separately. 

When the edge e is neither a bridge nor a looping edge, the deletion-contraction formula 
for the multivariate Tutte polynomial specializes to one for Pq: 

(3.7) P G (q,t) = P G ^ e (q,iM)+t e P G / e (q,$ e) ). 

In this case the numbers of connected components of G, G\e, G/e are all equal. Therefore, 
the same formula holds for Qq. It also holds once q is set to in the latter. 

We say that an edge is a regular edge if it is neither a bridge nor a looping edge. 

Proposition 3.8. Assume e is a regular edge of G. Then 

{W G } = L • [A* E ^ - (W G ^ e n W G/e )] - {W G/e } , 

{y G } = l • [A*^)- 1 - (y G ^ e n y G/e )] - {y G/e }. 

and the formula for {V G } is then the difference of these. 

Proof. The proof is entirely analogous to the one given for {Z G }, applied to the polyno- 
mials Q G and Q G \ q =o- □ 



The cases of bridges and looping edges are already dealt with in Lemma 3.7 we repeat 
them here for clarity: 

Proposition 3.9. If e is a looping edge of G, then 

{V G } = (T + l){V Gxe }. 

If e is a bridge, then 

{V G }-T{V G/e }. 
4. Edge splitting 

We now use the deletion-contraction formula of Theorem 13.11 to describe the effect of 
splitting an edge in a graph. 

Definition 4.1. Given a finite graph G and an edge e £ E{G), let °G denote the contrac- 
tion G/e; let 1 G = G; and more generally let k G be the graph obtained by replacing the 
edge e with a chain of k > 2 edges. 
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4.1. Splitting an edge. We have the following formula for the class of the hypersurface 
complement of the graph 2 G. Here and in the following, we denote by V(fi, . . . , fk) the 
zero locus {/i = • • • = fk = 0} of the ideal generated by (/i, . . . , fk). 

Theorem 4.2. For G a finite graph and 2 G the graph obtained by splitting an edge e in 
G, the class {Z2q} satisfies 

(4.1) {Z, G } = (T — 2){Z 1G } + (T - l){Zo G } + (T + 1) ({Z G ^ e } + {A G }) , 
where A e G = V(q + t ei Z G ^ e - qZ G/e ). 



The term {A G } appears to be difficult to evaluate geometrically. Equation 4.1 gives a 
relation between this term and the terms {Zmg}, which will allow us to obtain recursive 
formulas for these classes which are independent of {^4 G }. 

Proof. First observe that the effect of attaching an edge e to a vertex of a graph G is to 
multiply the polynomial Z G (q,t) by (t e + q) in the case of a non-looping edge, and by 
(t e + 1) in the case of a looping edge. This has the effect, in both cases, of multiplying the 
class {Z G } by the factor L — 1 as seen in Corollary 3.2 (5). 

Applying Theorem 3.1 to the graph 2 G requires handling the contraction, l G in this 
case, and the deletion, which is obtained from G\ e by adding a non-looping edge. In terms 
of equations, these are given respectively by the vanishing of Z G (q, t) and of Z Gxe (q, t^) ■ 
(q + t e ). The most interesting term is, of course, the class of the intersection of these two 
hypersurfaces, with ideal 

(4.2) (Z G ^ e (q,i^)(q + t e ),Z G (q,t)). 

This ideal defines a subscheme of A# E ( G *> +1 , while Zi G lives in A# B ( G ' +2 . Deletion- 
contraction applied to G gives 

Zg = Z G ^ e + t e Z G / e , 

therefore the zero locus V(Z G ^ e (q + t e ), Z G ) equals 

V(Z G ^ e (q + t e ),Z G ^ e +t e Z G/e ) = V{Z G ^ e ,t e Z G/e ) U V(q + t e ,Z G ^ e + t e Z G/e ) 

= V{Z G ^ e ,t e ) U V{Z G ^ e , Z G/e ) U V(q + t e , Z G ^ e - qZ G/e ). 

To apply inclusion-exclusion, we need the double and triple intersections of these compo- 
nents: 

V(Z G ^ e , t ei Z G ^ e , Z G / e ) = V(t e , Z G ^ e , Z G / e ) 

V(Z G ^, ei t e ,q + < e ,Z Gve - qZ G/e ) = V(q,t e ,Z G ^. e ) = V{q,t e ) 
V(Z G ^ e , Z G / e ,q + t e , Z G ^ e - qZ G / e ) = V(q + t e , Z G ^ e ,Z G / e ) 

and 

V{Zc^e, t e , Z G ^ e , Z G / e ,q + t e , Z Gve — qZ G / e ) — V(q, t ei Z G ^ e , Z G / e ) = V(q, t e ), 

where we used the fact that Z G ^ e and Z G / e are multiples of q. This implies that the triple 
intersection is in fact a double intersection, causing a useful cancellation at the level of 
Grothendieck classes: 

[V(Z G ^ e (q + t e ),Z G )} = [V(Z G ^ e ,t e )\ + [V(Z G ^ e , Z G/e )] + [V(q + t e , Z G ^ e - qZ G/e )\ 

- [V(t e , Z G ^ e , Z G/e )] - [V(q + t e ,Z G ^ e , Z G/e )}. 
All but one of the classes on the right-hand side have a clear interpretation: 

[V{Z G ^t e )\ = [Z Gve ], 
where we view Z Gxe as a hypersurface of A^ E ^; 

[v(z G ^ e , z G/e )] = l • [z G ^ e n z G/e ], 
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where the intersection is again viewed in A# E ( GS) and the factor of L is due to the free 
variable t e ; 

[V(t e ,Z G ^ e ,Z G/e )] = [V(q + t e ,Z G ^ e ,Z G/e )} = [z Gxe nz G/e ], 

still in A# E ( G )-. indeed, t e and q + t e may be used to eliminate t e in both cases, and both 
ideals define the same locus in A# E ( G > (with variables (q,i^)) after this projection. 
The remaining term is 

[V(q + t e , Z G ^ e - qZ G/e )} = [V(Z G ^ e - qZ G/e )}, 

where again we use q + t e — to eliminate t e , and view the right-hand-side as the class of 
a locus in A# E ( G \ with variables (q,i^). 

Let A G denote the locus determined by this ideal, as a subset of A# E ( G \ We have 
obtained that 

[V(Z G ^ e (q + t e ), Z G )] = (L - 2) ■ [Z G ^ e n Z G/e ] + [Z G ^ e ] + [A G ]. 
This can be equivalcntly stated in terms of classes of hypersurface complements as 

(4.3) {V(Z G ^ e (q + t e ), Z G )} = (L — 2) • {Z G ^ e n Z G/e } + {Z G ^ e } + {A G }. 
Indeed, we have 

{V(Z G ^ e (q + t e ),Z G )} = Ll £ l +1 - [V(Z G ^ e (q + t e ),Z G )} , 
while the complements of the other loci are taken in A^L Thus, 

{v(z G ^ e ( q + t e ),z G )} = i* E ^ +1 - ((l - 2) ■ [z Gxe n z G/e ] + [z G ^ e ] + [a g ]) 
= u^* E(G) - [z G ^e n z G/e ]) - 2(L# B < G ) - [z G _ n z G/e ])+ 2h# E ^ - {[z G ^ e ] + [A e G ]) 

with the stated result (4.3 1. 

Thus, we have obtained in this way an explicit calculation of the intersection term 
needed to apply Theorem 3.1 to the graph 2 G obtained by splitting an edge of G into two. 
We obtain 

(4.4) {Z2 G } = L((L-2){Z G ^ e nZ G/e } + {Z G ^ e } + {A G })-{Zi G } . 

We then apply again Theorem |3.1| to G to provide an alternative expression for the 
intersection term {Z G ^ e n Z G / e } and we obtain 

(4.5) L • {z Gxe n Z G/e } = {Z G/e } + {Z G } = {Zo G } + {Zi G }. 

□ 

The locus A e G determined by the ideal (q + t e , Z G ^ e — qZ G / e ) has an interpretation in 
terms of the combinatorics of the graph G. Let Z' denote the sum 

(4.6) e q k{A) n 

ACE aeA 

where the sum is restricted to the subgraphs not including e and connecting the endpoints 
of e; and let Z" denote the same expression, where the sum is restricted to the subgraphs 
not including e and not connecting the endpoints of e. Notice that there is a bijection 
between the monomial of Z G / e and the monomials of Z G ^ e ; in fact, 

(4.7) Zc^^Z' + Z" , Z G/e = Z' + ^-. 

Indeed, the graphs in Z" lose one connected component when e is contracted. 

Lemma 4.3. The locus A G may be described as V(q + t e , (l+t e )Z'), where the polynomial 
(1 + t e )Z' is the sum over all subgraphs of G (including e or not) which connect the 
endpoints of e in some way other than through e. 
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Proof. By (4.7) 



Z G ^ e - qZ G/e = (1 - q)Z' . 

Modulo q + t e this equals (1 + t e )Z' , and this term has the interpretation detailed in the 
statement. □ 

4.2. Multiple splittings. We can now formulate the result for multiple splittings of an 
edge e in a graph G, that is, for all the classes {Zmg}. 

Lemma 4.4. For all m > 1, the classes {Z m +i G } satisfy 

(4.8) {Z„» +1G } = (T - 2){Z mG } + (T — l){Z»-x G } + (T + l)T m_1 ({Z Gn J + {A G }) . 

Proof. If a polynomial is multiplied by t + g, where i is a new indeterminate, then the 
effect on the class {•} of the variety defined by that polynomial is to multiply it by 



L — 1 = T = [G m ] (cf. Corollary 3.2 (5)). Applying this observation to Z Gxe and A e G 
shows that 

{Z mGxe } + {A^ G } = T m ~\{Z G ^} + {A G }). 

where e denotes the last edge added in the splitting m G. Indeed, Z™ G ^ e is obtained from 
Z G ^ e by attaching a chain of to — 1 edges to G \ e: 

e 



and this has the effect of multiplying Z G ^ e by a term OIILi 1 + j where t\,... , i m _i are 
the variables corresponding to the edges in the chain. Therefore, {Zm G ^ e } = T m_1 {Z Gxe }- 



The effect on A e G is precisely the same. Indeed, recall that, by Lemma 4.3 the equation 
for A G in the hyperplane t e — —q is (1 — q)Z' ', where Z' is the sum corresponding to 
subgraphs of G \ e which connect the endpoints of e. From this description it is clear that 
the equation for A%. G in {t e = —q} is (1 — q)Z' J]™^ 1 (ti + q), and again it follows that 
{^4„ G } = T m ^ 1 {A G }. ' □ 



Combining Thcorcm |4.2| and Lemma |4~4| allows us to obtain a 'combinatorial' expression 
for the class {^4 G }. We then have the following recursive formula for the classes {Zm G }. 

Theorem 4.5. For all m > 0: 

(4.9) {Z m+3G } = (2T - 2){Z m+2G } - (T 2 - 3T + l){Z m+ i G } - T(T - l){Z mG }. 



Proof. Starting with (4.8), Theorem 4.2 allows us to express the term {Z G ^ e } + {A G } in 
terms of splittings: 

(4.10) (T + l)({Z G ^ e } + {A G }) = {Z. G } — (T — 2){Z 1G } — (T — l){Zo G }. 

Using the case to = 2 in Lemma |4.4| we then get 

{Zs G } = (T — 2){Z. G } + (T - 1){Z 1G } + T({Z 2G } — (T — 2){Z 1G } - (T — l){Zo G }), 
that is, 

{Zs G } = (2T - 2){Z 2G } - (T 2 - 3T + l){Zx G } - T(T - l){Zo G }. 



Then applying this formula to m+1 G rather than 1 G one obtains the recursion (4.9). □ 



One can then write a generating function for the classes {Zm G }, in a way similar to the 
corresponding result given in pQ for the graph hypersurfaces of Feynman graphs. 
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Theorem 4.6. The generating function of the classes {Zm G } is given by 

{Z» G } 



(4.11) "£{Z mG } 

m>0 

+ ( (T - 1 i 

Proof. Putting 
(4.12) 



JT-1). _ ( T _ i) 
e (T-l)-_, - 



T 



T + 1 

— (T — 2) 



e Ts - e- s 
' T + 1 

T-l)s _ e -s 
f 



T + 1 



{Z*g}- 



G e (s) := ]T{Z mG } 



m>0 



m! 



the recursion formula (4.9) translates into the differential equation 

(4.13) Gl'{s) = (2T - 2) G'^(s) - (T 2 - 3T + 1) G' e (s) - T(T - 1) G e (s), 
with solution 

(4.14) G e (s) = Ae- S + B e Ts + C e^~ 1)s . 

We can impose that this series begins with three undetermined coefficients, and solving 
for A, B, C gives 

{Z2 G } + {Z 1G } {Z2 G } + 3{Zi G } + 2{Zo G } 

^ = {Z ° g} + T f+1 

{Z 2G } + 3{Z lG } + 2{Zo G } 
T + l 

{Z2g} + {Zl G } 

T 

This then gives the form ( |4.11 ) for the generating function G e (s). 



B = -{Z 1G } - {Zo G } 
C = {Z 1G } + {Zo G } 



□ 



One can also write (4.111 in a form that involves explicitly the term A e G , with slightly 
simpler form of the coefficients, as follows. 



Corollary 4.7. The generating function (4.11| can be equivalently written as 

„Ts 



(4.15) Ei Z -o}-, = e- s 1 



m>0 



1 



T 



{Zo G } 



e (T+l)s _ e Ts 



e Ts -l 
T 

e Ts -l 



T 



{^g} 

({Z G ^ e } + {A G }) 



Proof. This follows directly from (4.111 and (4.10). 



□ 



4.3. Edge splitting for fixed q. The discussion is entirely parallel to the one given 
above for the class {Z G } with variable q. One has the analog of Theorem 4.2 in the case 
with fixed q, given by the following. 

Theorem 4.8. Let 2 G be the graph obtained by splitting an edge e in a graph G. Then 
the class {Z2 Gq } satisfies 

(4.16) {Z, G J = (T — 2){Zi G ,J + (T — l){Zo G ,J + (T + 1) ({Z G _, q } + {A G J) . 



Proof. As in the proof of Theorem |4.2| the main point is the computation of the class of 
V(Z G ^ e (q + t e ),Z G ). 
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Lemma 4.9. With q ^ fixed, the class of the locus V{Z G ^ e (q + t e ), Zq) is given by 
{V(Z G ^ e (q + t e ), Z G )} = (T - 1) • {Z G ^ q n Z G/e J + {Z G ^ q } + {A e G J, 

where A e G q C A* E< - G ^ 1 is the zero-locus of the polynomial Z G ^ e — qZ G / e = (1 — q)Z' , 
with fixed q, with Z' as in (4.6) (with fixed q). 



Proof. The argument here parallels closely the computation in the proof of Theorem 4.2 



In that computation (for variable q) there is a key cancellation of a class [V(g,i e )] due to 
inclusion-exclusion; for fixed q ^ this locus is empty. All other classes admit the same 
interpretation for fixed q as for variable q. 



In particular, the term A e G is the zero locus of (1 — q)Z' with Z' as in Lemma 4.3 (and 
q is now a fixed nonzero number). Note that this is when q = 1; in this case the last 
summand in the formula in Lemma [4~9| is 0. In general, arguing as in Lemma |4~3"1 one sees 
that Z' is given by the sum (4.6), over the range specified there. □ 

Implementing deletion-contraction gives then the same formula as in the case with 
variable q, so that one obtains the following. 

Corollary 4.10. Let 2 G be the graph obtained by splitting an edge e in a graph G. Then 
the class \Zi Gq \ satisfies 

(4.17) {ZzcJ = L ((L — 2){Z G ^ q n Z G/e J + {Z G ^ q } + {A G J) - ^ G , q ). 
Note again that the term {A e G } would be missing in the case q = 1. 



Next, use Proposition 3.4 to get a different expression for {Z G ^ e P\Z G / e }, and this gives 



a perfect parallel to Theorem |4.2| and completes the proof of Theorem |4.8| □ 



We can now pass to the case of multiple splittings, which again is analogous to the case 
with variable q. 

Theorem 4.11. Let m G be the graph obtained by multiple splitting on an edge e in G. 
Then, for all m > 0, the classes for fixed q satisfy 

{Z m+1G J = (T — 2){Z mGtq ) + (T - l){Z»-i G ,,} + (T + l)!™- 1 ({Z G ^ q } + {A G J) , 
which then gives the recursive relation 

(4.18) {Z m+3G J = (2T - 2){Z m+2G J - (T 2 - 3T + l){Z m+ i G , J - T(T - l){Z mG J, 

Proof. The argument is completely analogous to the case with variable q. The key step is 
the relation 

{A<L G J = T m - 1 {A G J, 
that one can see continues to hold in this case by the same argument used before. Then 



(4.18) is proved by using Theorem 4.8 to solve for the class ({Z G ^ eiq } + {A G }). □ 



The bottom line is that the same recursion holds for any fixed q as for the free q 
case. What will change are the seeds of this recursion, that is, the values of {Zm G q } for 
m = 0, 1, 2; these will naturally be different for fixed q. 

Nothing in the proof of the recursion excludes the case q = 1, and indeed {Zm G q=1 } = 
f#E(G)+m-i j g a so i u ti on of the recursion 

(2T - 2)T #B(G)+m+1 - (T 2 - 3T + i)t#- b (g)+™ _ j(j _ 1 -j T #-E(G)+m-i = t #e(g)+to+2^ 



With respect to the question of when the fibration condition (2.24) relating {Z G } to 
the {Z G q } holds, notice that, although the recursion is the same for all fixed q (all being 
the same as for the case of variable q) , one does not a priori know whether the seeds of 
the recursions are independent of q. 
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4.4. Edge splitting and the tangent cone. Again we consider first the operation of 
splitting one edge in two and then the case of multiple splittings. In the case of the tangent 
cone, we need to distinguish regular edges from bridges and looping edges. 

Lemma 4.12. If the edge e of G is either a bridge or a looping edge and m G denotes the 
graph obtained by iterated splitting of the edge e, then for m > 1 the class of the tangent 
cone complement satisfies {Vm+i G } = T m {V G }. 

Proof. If e is a bridge, then splitting e amounts to inserting a new bridge; by the deletion- 
contraction formulas for bridges (Proposition 3.9 ) 

{V2 G } = T ■ {V G }. 

Thus, in the case of a bridge one also already sees that splitting the edge multiple times 
just has the effect of multiplying {V G } by a power of T. 

If e is a looping edge, adding a loop to a graph G' multiplies the class {Vc} by T + 1 
as seen in Lemma 3.7 Attaching a split loop amounts to multiplying Zq> (and hence P G >) 
by (q + t e + tf), where t e and tf are the variables corresponding to the two new edges; 
this is simply because q 2 + t e q + tjq is the ^-polynomial for a 2-banana. Here G' is the 
graph obtained from G by removing the loop. From 

P2 G = P G , (q + t e +t f ) 

we see that Pi G ^ implies P G i ^ and t e ^ — (q + tf). 

Therefore, for any point (q, t) for which Pqi ^ we have an A 1 worth of choices for tt 
and an A 1 \ A worth of choices for t e . So we can conclude that 

{W2 G } = L(L - 1) ■ {W G >} = T • {W G }. 

The same analysis goes through after setting q to 0, and hence 

{y* G } = l(l - 1) • {y G >} = t ■ {y G }. 

The conclusion is that 

{V 2G } = T • {V G }. 

Thus, splitting a loop once has again the effect of multiplying the corresponding {V G } by 
T, as for bridges. □ 

We then check the more interesting case of regular edges. 

Theorem 4.13. Let e be a regular edge of G, and let 2 G be the graph obtained by intro- 
ducing a 2-valent vertex in e, thereby splitting it. With other notation as above, 

{V2 G } = (T — 2){V G } + (T — l){V G/e } + (T + l)({V G ^ e } + MQ Gxe - qQ G /e)})- 

Proof. Let e be a regular edge, and call e (again) and / the two edges created in the 
process. The polynomial for Wi G is found easily by applying deletion-contraction: 

Q 2 G = Qc\e ■ (q + te) + tf ■ Q G . 

Indeed, deleting / leaves the graph G\e with a dangling edge e attached, and attaching an 
edge e to a vertex has the effect of multiplying the corresponding polynomial by (q + t e ): 
contracting / gives G back. By Proposition |3.8| therefore we have 



(4.19) {W2 G } = L • [A# E ^ +1 - W n ] - {W G }, 

where Wn is the intersection in A^ E ^ +1 of the hypersurface W G , with equation Q G = 0, 
and the hypersurface with equation Q G ^ e (q + t e ) = 0. 
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Set theoretically, the locus Wn is given by 

Wn = V(Q G ^ e {q + t e ),Q G ) = V(Q G ^ e (q + t e ), Q Gxe + t e Q G/e ) 

= V{Q G ^ e , t e Q G/e ) U V(q + t e , Q Gxe + t e Q G/e ) 

= V(Q G ^ e , Q G/e ) U V(t e ,Q G ^ e ) U V(q + t e ,Qa^e - qQc/e), 

where all ideals are viewed in A# E ( G ^ +1 , with coordinates q and t a , a € E(G). Inclusion- 
exclusion then gives [Wn] as the sum of the three loci indicated here, minus the sum of the 
three pairwise intersections, plus the triple intersection. The three pairwise intersections 
are 

V(t e , QG X e, Qg/c) 
V(q + t e ,Q G ^ e ,Q G/e ) 
V(q,t e ,Q G ^ e ) 

and the triple intersection is 

V(q,t e ,Q G ^ e ,Q G / e ) . 
This is very similar to the situation we have seen for the case of Z G . We find 

[V{Q G ^e, Q G /e)} = L • [W G ^ e H W G/e ] 

as G \ e and G/e have E(G) \ {e} as index set, while 

[V(t e , Qgxc QG/e)] = [V(q + t e , Q G \e, Qc/e)] = [W G ve D W G /e] ■ 

In both cases, the first equation eliminates t e , doing nothing to the rest since the rest does 
not depend on t e . We also have 

[V(t e ,Q G ^ e )] = [W G ^ e ] , [V(q,t e ,QG^e)} = [y G ^e\, 

and 

[V(q,t e ,Q G ^e,Q G /e)] = [y G ^eny G/e ]. 

That leaves us with 

V(q + t e , Q G ^ e - qQ G/e ) C K# E ^+\ 

or simply 

V(Q G ^-qQ G/e )Qk# E ^ G \ 
since the only effect of the first equation is to eliminate t e . 
Applying then inclusion-exclusion we obtain 

[W n ] = (L — 2)[W G ve n W G/e ] + [W G ^ e ] + [V(Q G ^ e ~qQ G/e )} - [y G ^ e ] + [J^Gve H ^/J . 

Passing to the classes of the complements and implementing the resulting expression for 
2 G then gives 

{W2 G } = (L — 2)U^* E[G) - [w G ^ e n W G/e ]) + L{W G ^ e } 

+ L{V(Q G ^ e - q Q G / e )} - HyG^e]} + H^* E{G) ~ X - t^CKe D Jfe/el) - {W G }. 

Using Proposition 3.8 and computing the difference {V G } = {W G } — {y G } one gets 

{V2 G } = {m G } - {y* G } 

= (l - 3){w G } + (l — 2){w G/e } + L({w Gxe } - {y G ^ e }) + {y G } 

+ {y G /e} + HV(Q G ^e - q Q G /e)} ~ (L - 2){^ G } — (L — l){y G/e } - L{}^ e } 

= (l - 3)({w G } - {y G }) + (l - 2)({w G ^ e } - {y G ^}) 

+ L({W G -. e } - {y G ^ e }) + H{V(Q G ^ e - q Q G /e)} - {y G ^e}) 

= (L - 3){V G } + (L — 2){V G/e } + L{V G ^ e } + h{V(Q G ^ e - qQ G/e )}. 

□ 
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As in the case of Z G and Z Gq , the polynomial Qg xe — qQ G / e can be given a combi- 
natorial interpretation as 

QGve — q Qc/e — £ q k ^- k ^\[t a . 

A C E \ {e} forest connecting endpoints of e (iGA 

We then can proceed as in the case of Z G and Z G ,q and obtain the following recursions. 
Proposition 4.14. With notation as above, and assuming e is not a looping edge, 

{y™+z G } = (2T - i){y m+ 2 G } - t(t - 2){y m+ i G } - r 2 {y mG } 

for all m> 0. The generating function 

m>0 

/or i/ie classes {3^g} *s iften o/ £/ie form 

G*(s) = Ae- S +Bse Ts +Ce Ts , 
with the constants A, B, and C satisfying 

,{^ G } + {3M , {y»G} + 2{yi G } + {y° G } 



A = {y 0G }-2- 



b = -{yi G } - {y° G } + 

C = 2 



T + l (T + l) 2 

{y> G } + 2{y 1G } + {y 0G } 



T + l 

{yi G } + {yo G } {y* G } + 2{y 1G } + {y 0G } 



T + l (T + l) 2 

Similarly, 

{Vm + a G } = (2T - 2){V™ + 2 G } - (T 2 - 3T + l){V m+1G } - T(T - l){V mG }, 
for all m > 0, with generating function 

G v e (s) : =]T{y mG }: 



s 

ml 

m>0 



given by 

Gg (s) = Ae~ s + B e Ts + C e^ T_1 ' s 
with the terms A, B and C satisfying 

{V 2G } + {V 1G } {V 2G } + 3{V 1G } + 2{Vo G } 



A = {Vo G } 



B = -{V 1G } - {Vo G } + 
C = {V 1G } + {Vo G } 



T T + l 

{V 2G } + 3{Vi G } + 2{Vo G } 



T + l 

{v^ G } + {y 1G } 

T 



The argument is essentially analogous to the cases of Z G and Z Gtq analyzed before and 
we do not reproduce it explicitly here. 

The generating function for y™. G is 'dual' to the one for X™ G given in [T]: the effect of 
splitting edges on the class {3^ m G } is analogous to the effect of multiplying edges on the 
class {X™ G }. 
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5. Polygons and linked polygons 

We now focus on a particularly simple class of graphs for which we can compute ev- 
erything explicitly. These will be polygons and graphs constructed out of chains of linked 
polygons. We will later focus especially on this class of graphs to provide an explicit ex- 
ample of how to apply these motivic techniques to analyze (virtual) phase transitions in 
the corresponding Potts models. 

We start by using the formulae for edge splittings obtained in the previous section to 
compute explicitly the classes {Z G } for polygon graphs. 

Proposition 5.1. Let m G be an (m + 1) -sided polygon. Then the classes {Zm G } are given 
explicitly by the formula 

("TP — 1 Y m — (— i V™ 



(5.1) {Z mG } = T m+2 +T(T-l)(T m -(T-l) m ) + (T-l)- 

Proof. Let us check directly the initial cases that are needed to use the recursive formula 
for edge splittings. The first graphs are °G a single loop, 1 G a 2-banana (two vertices with 
two parallel edges between them), 2 G a triangle. The equations Z G (q,t) = are of the 
form 

°G : q + qt = (in A 1+1 ) 

l G : q 2 + (h + t 2 + t^q = (in A 2+1 ) 

2 G : q 3 + (h + t 2 + t 3 )q 2 + (^ 2 + ht 3 + t 2 t 3 + ht 2 t 3 )q = (in A 3+1 ) 

The corresponding classes {Z G } can then be computed directly in these cases by applying 
the basic facts listed in Corollary |3 . 2 1 and Theorem |4.2| One obtains 

{Zo G } = T 2 

{Zi G } = T 3 +T 2 - 1 

{Z2 G } = T 4 + 2T 3 - 2T 2 - 2T + 2. 



The expression (5.1 1 is then ml times the coefficient of s m in the expansion of the right- 
hand side of the generating function (4.11 1, with these initial conditions. □ 

5.1. Polygons at fixed q. We will also need in the following the classes of the polygon 
graphs for fixed value of q. These are obtained as follows. 

Proposition 5.2. Let m G be an (m + l)-sided polygon. Then the classes {Zm G q } for 
fixed q ^ 0,1 are given explicitly by the formula 

(T — 1 ) m — (— 1 ) m 
(5.2) {Z mG: J = T m+1 +T(T m -(T-l) m ) + ^ ' 



T 



The polygon graphs m G satisfy the fibration condition (2.24). 



Proof. As above, the seeds of the recursion are a single loop, a 2-banana, and a triangle. 
The single loop has class T from Corollary |3.5| 

The 2-banana hypersurface has equation qt\t 2 + qt\ + qt 2 + q 2 = 0. Factoring out a q 
(assumed to be nonzero to begin with) this is equivalent to 

(ti + l)(t 2 + l) = l-g. 
We are assuming that q ^ 1, so up to a variable change this equation is 

u\u 2 = r^0. 

This forces u\ ^ 0, and determines u 2 once u\ is fixed; thus the class of this locus in the 
Grothendieck group is L — 1; its complement in A 2 has class 

L 2 - L + 1 = T 2 + T + 1. 
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This is independent of q ^ 0,1. 

The triangle hypersurface has equation 

q(tit 2 t 3 + ht 2 + htg + t 2 t 3 + q{ti +t 2 + t 3 ) + q 2 ) = 0. 

Changing variables: m = ti + 1, r = q — 1, and keeping in mind g 7^ 0, this is equivalent to 

(5.3) uiu 2 u 3 + (ui + u 2 + u 3 )r = r(l - r). 

Solving for u 3 gives 



u 3 



r(l — r — u\ — u 2 ) 



uiu 2 + r 

Thus, the variety contains an open subvariety isomorphic to 

A 2 \ V(uiu 2 + r) ; 

this locus has class L 2 — L + 1, as we just computed above, since r 7^ (as q ^ 1). If 
Uiu 2 — —r, then (5.3 1 is equivalent to 

ui + u 2 = 1 - r, 

implying easily that (ui,u 2 ) equals (l,r) or (r, 1), while u 3 is free in this case. That is, 
the complement in (5.3) of the open subvariety determined above consists of the loci 

{(T-r,^)}^ 1 , {(-r,l,u 3 )}=A 1 . 



These lines are distinct, since (as q 7^ 0). The conclusion is that the class of (5.3) 

equals 

L 2 -L+1 + 2L = L 2 +L+1 , 

and hence its complement in A 3 has class 

L 3 - L 2 - L - 1 = T 3 + 2T 2 - 2 . 

Again this is independent of q. 

These cases are compatible with the fibration condition (2.24). For instance, in the 
triangle case, we get 

(T - 1)(T 3 + 2T 2 - 2) + T 3 = T 4 + 2T 3 - 2T 2 - 2T + 2 

in agreement with the class for a triangle in the free q case, used in Proposition |5.1| 
The recursion then gives the classes (5.2), 



{Zm G<q } 



T m+2 +T(T - l)(T m - (T - l) m ) + (T - 1) ( T - 1 ) m -(- 1 )'" _jm+i 



T m+i +T(T m -(T-l) m ) + 



T - 1 

(T - l) m - (-1) 7 
T 



One can verify explicitly the compatibility of (|5.2|) and (|5.ip with the fibration condition 

(T - l) m - (-1)" 1 



(2.24): the polynomial 

{Z™ G } = T m+2 + T(T - l)(T m - (T - l) m ) + (T - 1) 



T 



factors exactly as predicted by (2.24) in terms of the classes (5.2). 



□ 
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5.2. Polygons and the class of the tangent cone. The explicit formula for the classes 
{V G } of the complement of the tangent cone of the Potts model hypersurface for polygons 
is obtained as follows. 

Proposition 5.3. Let m G be the polygon with m + 1 edges. Then 

{V mG } = (T - l)(-l) m + 2T m+2 - (T + 1)(T - l) m+1 . 
Proof. For the seed of the recursion we have in this case 
°G : q = (inA 1+1 ) 

,2 , (• , . \ n /• a 2+1 



2 G 



T I ' T 



<r + (*i+ta)? = (inA 2+i ) 

q 3 + (i x + i a + t 3 )q 2 + (tit 2 + t x t 3 + t 2 t 3 )q = (in A 3+1 ) 

for which we then get 

{Vo G } = L 2 - L 1 =T(T + 1) 

{Vi G } = L 3 - 2L 2 + L 1 = T 2 (T + 1) 

{V2 G } = T 4 + 2T 3 - T = T(T + 1)(T 2 + T - 1). 

The first two expressions are immediate; for the third, note that V2 G is a cone over a 
nonsingular quadric in P 3 ; the computation is then straightforward. We then have 

GX{s) = e- s ( (l + T(T + 1) + ^^T 2 (T + 1) 

e (T+i )s _ e Ts _ \ . 2J : 

that is 

G^(s) = (T — l)e- s + 2T 2 e Ts - (T 2 - l)e (T - 1)s . 
Reading off the coefficient of s m /ml then gives the result. □ 

5.3. Chains of linked polygons. A class of graphs that have been extensively studied 
from the point of view of the statistical mechanics of Potts models is the case of chains of 
linked polygons. These are graphs consisting of N equal polygons rn G, each with m + 1 
edges, attached to one another by chains of k > edges. The case k = corresponds to 
polygons joined at vertices. The corresponding Potts models were studied, from the point 
of view of the properties of ground state entropy, in [52] : see also the references therein 
for several other results on this class of Potts models. 




Definition 5.4. Let ( m > k )Q N denote the graph obtain by joining N polygons, each with 
m + 1 sides, with every two nearby polygons connected by a chain of k > edges. 

This family of graphs has three parameters m, k, and N, each of which can be inde- 
pendently sent to oo to create an infinite graph. Thus, we will use this family as our main 
example on which to analyze the topological complexity of the corresponding set of virtual 
phase transitions. 

The classes {Zq^} can be computed explicitly for this family from the result for poly- 
gons in Proposition |5.2| and the other basic properties of Corollary |3.5| 
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Proposition 5.5. For^ m ^G N as in Definition\5.A the classes {Z Gq \ with fixed q 7^ 0, 1 
are given by 

(nr _ i\m _ ( _-| \m\ N 
pm+1 _|_ fjp/fjpm /rjp l) m ) -4- — - I ^h(N— 1) 



(5.4) ( T m+i + T(T m - (T - l) m ) + 



Proof. By (3) and (4) of Corollary 3.5 we have 



{z (m , k)GN j = {z,^ q } N Tr fe ( JV_1 ). 



The result then follows from (5.2). □ 



In this case, since the graphs involved satisfy the fibration condition (2.241, one can 
also obtain an explicit formula for the classes {Z( m ,k) G N} with variable q, from (3.3) of 
Corollary |3.3| 

6. Multiple edge formula 

The operation of doubling edges is dual to splitting edges. However, while in the case 
of the graph hypersurfaces of Feynman graphs analyzed in [I] the corresponding operation 
on the hypersurface complement classes in the Grothendieck ring is very simple, this is 
not the case when one considers the Potts model hypersurface. 

In fact, the combinatorial deletion-contraction formula for the multivariate Tutte poly- 
nomial shows that the polynomial for the graph obtained from G by doubling an edge e 
is of the form 

(6.1) Z G ^ e + {t e + t f +t e t f )Z G/e , 

and it is the presence of the extra term t e tf here that complicates the matter. 

6.1. Edge doubling. We derive here a formula for the class {Z G } under the operation of 
doubling an edge, and then we obtain a recursive formula for the iteration of this operation. 

Theorem 6.1. Let G' be the graph obtained by doubling the edge e in a graph G. Then 

(6.2) {Z G ,} = T.{Z G } + (T+1)-{B G }, 
where B G is the locus of zeros of Zo^ e — Z G / e . 

Proof. It is convenient to change variables, letting u e = 1 + t e , Uf = 1 + tf. Then the 



class (6.1 1 for the double edged graph is given by 

(6.3) Z G ^ e + (u e Uf - l)Z G/e ^ 0. 

If Z G / e = 0, then necessarily Z G ^ e ^ 0; u e and Uf are free, so this accounts for a class 

(T + i) 2 -[z G/e \(z Gve nz G/e )]. 

If Z G / e ^ 0, then the condition amounts to 

, 1 Z G ^ e 
u e u f £1- — . 

Z G/e 

This in turn leads to to two possibilities: 

• Either = 1, and then u e Uf =^ 0; this accounts for L 2 — 2L + 1 = T 2 ; 

• Or Z z ^° 7^ 1, and then u e Uf ^ c for some c ^ 0. For c ^ 0, u e Uf = c necessarily 

gives Ut 7^ 0, u e — c/uf. this accounts for L — 1. Thus the class of u e Uf ^ c is 
L 2 - L+ 1 = T 2 +T + 1. 
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In total, the class of the complement for the double-edged graph is 

(T + 1) 2 • [Za/e \ [z G ^ e n z G/e )} + T 2 [(Ai £ i \ z G/e ) n v(z G ^ e - z G/e )] 

+ (T 2 + T + l)[(Al £ l n Z G/e ) \ V(Z G ^ e ~ Z G/e )]. 
This expression can be further simplified in the following way: 

(T + l) 2 • [Z G/e \ (Z G ^ e n Z G/e )} + (T 2 + T + 1)[A* E ^ \ Z G/e ] 

- (T + 1)[(A* E ^ \ z G/e ) n V(Z G ^ e - Z G/e )] 

= (T 2 + 2T + 1) ■ [Z G/e ] - (T + l) 2 [Z G ^ e n Z G/e ] + (T 2 + T + 1)[A* E ^ \ Z G/e ] 

- (T + 1)[(A* E ^ v Z G/e ) n V(Z G ^ e - Z G/e )] 

= T • [Z G/e ] - (T + l) 2 [Z G ^ e n Z G/e ] + (T 2 + T + 1)L #B ( G ) 

- (T + 1)[(A# B < G ) v Z G/e ) n V(Z G ^ e - Z G/e )} 

= (T + 1) 2 ([L# £ ( G ) - [Z G _ n Z G/e ]) - T([h* E ^ - [Z G/e ]) 

- (T + 1)[(A#^ G > \ z G/e ) n V(Z G ^ e - Z GJe )]. 



We then use Theorem 3.1 to express {Z G ^ e PI Z G / e }, obtaining 

(T + 1)({Z G } + {Z G/e }) - T({Z G/e }) - (T + 1)[(A# B < G > \ Z G/e ) n ^(Z G ^ e - Z G/e )] 
and hence simply 

(T + 1){Z G } + {Z G/e } - (T + 1)[(A#^ G > \ Z G/e ) n V(Z Gve - Z G/e )]. 

Now 

( A #E(G) x Zc/e) n F(Zg ^ £ _ Zc/e) = y(Zc ^ e _ ZG/e) x (Zgx£ n Zc/e)j 

and again using Theorem |3.1| we get that the class of the complement for the double-edged 
graph is 

T • {Z G } - (T + 1) • {V(Z G ^ e - Z G/e )}. 

□ 

As in the case of the term A G in the formula for edge splitting, the term B G here also 
has an interpretation in terms of the combinatorics of the graph. 

Lemma 6.2. The locus B G of zeros of Z G ^ e — Z G / e is equivalently the locus of zeros of 

a polynomial (q — 1)Z , where Z" is the sum of the monomials IlceA^ 6 corresponding to 
subgraphs A of G/e that acquire an additional connected component when they are viewed 
in G \ e, and Z = Z" / q. 



Proof. Recall from (4.7 1 that 

Z" 

Z G ^ e = Z + Z , Z G / e = Z H . 

q 

and hence 

Z G ^.e - Z G / e = (q - l)Z , 

where z" = Z" jq is indeed the sum of the standard monomials over the subgraphs of G/e 
which acquire an additional connected component when they are viewed in G \ e. □ 
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Thus, the description of B e G is somewhat complementary to that of A e G : both are (q — 1) 
times a sum of terms having to do with subgraphs of G/e. For A, one looks at graphs 
for which the number of connected components is the same when the subgraph is viewed 
in G \ e; for B, one looks at the graphs for which the number of connected components 
increases. 

Example 6.3. When e is a looping edge, then B G — 0: indeed, G \ e = G/e in this case. 
Thus {B G } — 0, and the formulas simplifies to {Zc>} = T • {Z G }, which recovers the case 
(5) of Corollary 3.2 namely attaching a new looping edge to G. 

6.2. Multiple edge formulas. We now consider the case where the operation of doubling 
an edge is iterated, that is, where an edge e in a graph is replaced by m parallel edges, 
between the same endpoints. This can be seen also as replacing an edge e with the m-th 
banana graph (see [3]). 

Theorem 6.4. Let G^ m > denote the graph obtained by adding m edges parallel to e in G. 
(So G = G(°\ G' = GW.J Then, for m > 0, the classes {Z G(m) } satisfy 

(6.4) {Z G(m+2) } = (2T + l){Z G(m+1) } - T(T + l){Z G(m >}. 

Proof. The key case to consider is that in which we triple a given edge of G: let G' denote 
(as in Proposition 6.1 ) the graph obtained from G by doubling e, and let G" be the graph 



obtained from G' by doubling e again. Applying Theorem |6 . 1 1 yields 

{Z G „} = T-{Z G ,} + (T + 1).{B G ,}. 



Thus, we have to understand {B G ,}. According to Lemma 6.2 this hypersurface has 
equation (q— l)Z , where Z collects monomial according to the subgraphs of G'/e which 
acquire a component when viewed in G' \ e. The new edge in G' parallel to e cannot 
be part of any such subgraphs, since it does join the endpoints of e, so it prevents a new 
component from forming as we remove e. Therefore, the Z for G actually equals on the 
nose the Z for G; the only difference between B G and B G , is that the latter is contained 
in a space of dimension one higher, and it may be described as a cylinder on B G . 
Thus, we have 

{Z G „} = T-{Z G ,} + (T + l)-{B G ,} 
and {W G ,} = (T + 1){W G }. By Theorem |6~Tj 

(T + 1){W G } = {Z G ,}-T{Z G }. 

Thus, we obtain 

{Z G „} = T{Z G ,} + (T + 1)({Z G ,} - T{Z G }) 
= (2T + 1){Z G ,}-T(T + 1){Z G }. 

The stated formula follows by applying this formula to G'™-' instead of G. □ 

We can form a generating function for the classes of graphs with multiple edges. 
Theorem 6.5. The generating function of the classes {Z G ( m }} is given by 

(6.5) Yl i Z G^>}^ = ((T + 1){Z G } - {Z G ,}) e T * + ({Z G ,} - T{Z G }) e^ . 
^— ' m! 

m>0 



Proof. The recurrence relation ( 6.4 ) of Thcorem|6.4|translates into the differential equation 



(6.6) g"(s) = (2T + l)g'(s) - T(T + l)g(s) 

for the generating function 

(6-7) 9(s) = Y,{ z g^Y-i> 
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solving which reveals that 

Y^{Z G{m) } S —=Ae^+B^s. 
L — ' m! 



Solving for the constants A, B in terms of {Z G } and {Z G i} gives (6.5). □ 

Example 6.6. The m-th banana graph is a graph with two vertices and m parallel edges 
between them. To compute the class for bananas, we can start with G a single non-looping 
edge, for which {Z G } = T 2 , and G' a 2-banana, for which {Z G ,} = T 3 + T 2 - 1. Then 
from (6.5) we have 

E{ Z G (.")}^= e Ts + (T 2 -l)e (T+1)s , 
^— ' mi 

m>0 

from which we obtain 

(6.8) {Z Gim) } = T m + (T — 1)(T + l) m+1 

for the class of the m-th banana graph. A generating function for the analogous class for 
graph hypersurfaces X G is given in pQ. 

6.3. Multiple edges for fixed q. Again the argument for variable q carries over almost 
identically to cover the case with fixed q ^ 0,1. We obtain the following results. 

Proposition 6.7. Let G' be the graph obtained by doubling the edge e in a graph G. Then 

{Z G , tq } = T.{Z G , q } + (T + l).{W G>q }, 

where the locus W Gq C A*^ 6 )" 1 , for q^Q,l, is given by the vanishing of the polynomial 
Z" adding monomials over the subgraphs of G/e which acquire an additional connected 
component when they are viewed in G \ e. 

Corollary 6.8. Let G*- 1 ™-' denote the graph obtained by adding m edges parallel to e in G. 
(So G = G(°) , G" = GW Then for m > 

{Z G(m+2) , J = (2T + l){Z Glm+1) J - T(T + l){Z Gim)>q ). 

The general solution of this recursion also matches the one for free-g: 

m 

E i Z G^, J^j = « T + 1 ^{ Zg i} - i Z G', q }) e Ts + ({Za>, q } - T{^ G ,J) e (T+1)s . 

m>0 

Example 6.9. Consider the case of the banana graphs. The seeds for bananas are a single 
non-looping edge, with class T and the 2-banana, with class T 2 +T + l, as computed above. 
Plugging into the last formula, we get the generating function for the classes of Potts model 
complements of banana graphs for fixed q: 



' H ml 
m>0 

extracting the term of degree m gives the very simple class for the (m + l)-banana: 
(6.9) (T + l) rn+1 - T m 



in agreement with the fibration condition (2.24) for (6.8), and in particular independent 
of q ^ 0,1. 
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6.4. Chains of linked banana graphs. By analogy to the example of the chains of 
linked polygon graphs considered in |5.3[ we consider here a similar family but with the 
polygons replaced by banana graphs. 




Definition 6.10. The graphs k G^ m ^' N are obtained by connecting N banana graphs G^ m ' , 
each with m parallel edges, each connected to the next by a chain of k > edges ( connected 
by joining vertices in the case k = 0). 

We can compute the classes {Z Gtq } for this family of graphs using the explicit formula 



(6.9) for the banana graphs. 

Proposition 6.11. Let k G^ T 
classes are given by 



.N 



be the graphs of Definition 6.10 Then the corresponding 



(6.10) 



{Z fcG(m >,„ J = ((T + l) m+1 - jrnjN T k(N-l)_ 



Proof. The result follows immediately by applying (3) and (4) of Corollary 3.5 to the 



explicit formula (6.9) in Example 6.9 



□ 



6.5. Polygon chains. A class of graphs that can be obtained by alternating edge splitting 
and edge doubling operations in different orders are the chains of polygons of various sizes 
joined along edges. 



In the case of the graph hypersurfaces of Feynman graphs, explicit formulae for the 
Grothendieck classes of the hypersurface complements for this type of graphs were obtained 
in pp. However, in the case of the Potts model hypersurfaces, the recursion relation 
becomes more complicated: if one combines the formula for doubling an edge with the 
formula for splitting it, the resulting class becomes of the form 

(T 2 - T + 1){Z G } + (T - l)T{Z G/e } + (T + \){C C G ) + (T - 2)(T + l){D e G }, 

where the term C G is obtained combinatorially from the subgraphs of G that connect the 
endpoints of e, and D G is obtained from paths in G\ e which do not connect the endpoints 
of e. There is then no obvious recursive procedure of the type used for either the splitting 
or the doubling alone, which takes care of eliminating both of these additional terms. This 
remains an interesting case of graphs to investigate. 
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7. ESTIMATING TOPOLOGICAL COMPLEXITY OF VIRTUAL PHASE TRANSITIONS 

We show here, in the concrete example of the chains of linked polygon graphs ( m > k )Q N 
and in the simpler example of the banana graphs, how one can use our calculations of 
classes in the Grothendieck ring to estimate how the topological complexity of the set of 
virtual phase transition grows as the graphs grow in size within the given family. 

Good indicators of topological complexity are homologies and cohomologies and the 
associated Euler characteristics. In the case of the real algebraic varieties Za tq (M,), which 
can in general be singular and non-compact, therefore, we can take as a good indicator 



the Euler characteristic with compact support Xc(^G,q(^-))- As we discussed in ^ 2.8 this 
is the unique invariant of real algebraic varieties that is both a topological invariant and 
a motivic invariant. Using the fact that it factors through the Grothendieck ring T^oC^r)' 
we obtain the following results. 



Proposition 7.1. Let ( m > k )Q N ft e the chain of linked polygons graphs of Definition 5.4 
Then the Euler characteristic with compact support Xc(Z( m ,k) G N of the set of virtual 

phase transitions Z( m ,k) G N q (M) of the Potts model is given by 

^_^mN+kN-k ff_-^N _ 2 k N-k~N ^3™+! _|_ ^ _ 2 m +'^ N ^j 



Proof. As we have seen in Example 2.8 for real varieties Xc(T) = Xc( < G m (M)) = —2. Using 



the fact that \c is a rm g homomorphism Xc '■ -KoC^r) — > Z we then obtain from (5.4) of 
Proposition |5.5[ 

(— \y nN + kN - k 2 kN - k ~ N ^3 m + 1 + 1 — 2" l+3 ) JV 
We then use additivity again to obtain 

Xc{Z {m , k)GN , ,(M)) = X c{h* E( ^ k)GN) ) - Xc({iW)G", 9 }), 
where 4E{( m ^G N ) = N(m+ 1) +k(N- 1) and Xci^ 1 ) = Xc(R) = -1, so that we obtain 



(7.1). □ 



We obtain a similar result for the class of graphs k Q( m )' N of Definition 6.10 obtained 
by chains of linked banana graphs. 



Proposition 7.2. Let k Q( m )^ N he the graphs of Definition 6.10 Then the Euler charac- 
teristic with compact support Xc(Zk G ( m ),N (R)) of the set of virtual phase transitions in the 
model is given by 

(7 2) [_]_ynN+kN+N-k A _ 2*(W-1) (^m + \} N ^ 



Proof. The argument is exactly as in the previous case, using the expression (6.10) for the 



classes in the Grothendieck ring. □ 

7.1. Decision complexity and topological complexity. There is, in fact, a more 
technical sense, from the point of view of complexity theory, according to which the Euler 
characteristic Xc(Z G ^ q (M.)) really gives an estimate of the complexity of the real algebraic 
variety Z Giq (M.) of virtual phase transitions of the Potts model over the graph G. We 
describe it briefly following the survey [10] . 

An algebraic circuit over R is an acyclic directed graph where each node has in-degree 
either or 1, or 2. The nodes of in-degree are the input nodes and they are labelled 
by real variables; the nodes of in-degree 1 are either output nodes (out-degree equal to 
zero) or sign nodes: these are nodes that, to an input x assign output 1 if x > and zero 
otherwise; and the nodes of in-degree 2 are labelled by an operation +, — , x, or / and are 
called arithmetic nodes. The size cr(C) of an algebraic circuit C is the number of nodes. To 
an algebraic circuit C with n input nodes and m output nodes one can associate a function 
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tfc ■ K ra — ► R m > the function computed by the circuit. A decision circuit is an algebraic 
circuit with only one output node returning values of or 1. To each decision circuit 
one can associate a real (semi) algebraic set Sc — {x G M. n | (fc{x) — 1} and, conversely, 
it is known that each (semi)algebraic set in M. n is realized by some decision circuit. The 
decision complexity of a real (semi)algebraic set S is defined as 

(7.3) C(S) = mm{a(C) | S c = S}. 

It is known by [27] that there is a lower bound on the decision complexity of a (semi) algebraic 
set S of the form 

(7-4) C(S)> ^(log 3 x c (5)-n-4), 

in terms of the Euler characteristic with compact support. A similar, more refined estimate 
exists in terms of the sum of the Borel-Moore Betti numbers. 

Instead of considering the real algebraic varieties of virtual phase transitions Zq , 9 (IR), 
one can look at the actual physical phase transitions. For a finite graph G, this means 
considering, in the antiferromagnetic case, the semialgebraic set given by the intersection 
of Z Gi9 (E) with the semialgebraic set S = {t e A# E( - G ^> \ - 1 < t e < 0}. 
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